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SECTION 1 


INTRODUCTION AND SUMMARY 


1. 1 BACKGROUND 

The pressing need to better survey and manage the earth's resources and environment has 
prompted man to explore the possibilities of remote sensing from space. Early efforts 
began with space photographs from the Gemini and Apollo programs and continued with multi- 
spectral data from Landsat 1 and 2 spacecraft. Landsat D is currently planned as the ne.xt 
majoi step for the Earth Resources Program. 

Landsat 1, launched in 1972, marked the start of NASA's Earth Resources satellite program. 

This successful sixacecraft was followed two and a half years later with Landsat 2, an iden- 
tical spacecraft. The overwhelming success of these two Landsats, demonstrated through 
hundreds of exjjerimental programs, has motivated NASA to continue to improve the Earth 
Resources satellite program. The third satellite, Landsat C, has teen procured and is scheduled 
for launch in late 1977. This third satellite wdll carry a modified Multispectrul Scanner and 
will utilize an improved digital ground system. NASA is now planning for the next step, Landsat 
D, which will provide several major advances. Landsat D will incorporate the ITiematic 
Mapper (TM) as a new' sensor, it will utilize the Multi-mission Modular Spacecraft (MMS), 
it will make use of the Tracking and Data Relay Satellite System (TDRSS) and it will employ 
a new more advanced ground system. Each of these represent significant improvements in 
the sUite-of-the-art. This study is one of several which address various aspects of the planned 
I.andsat D system. 

As the Earth Resources Program has matured through the landsat spacecraft it has begun 
the transition from an experimental research activity to a sound demonstration of proven 
utility. This imptirtant transition will be completed with the landsat D system which in- 
corporates several key improvements over the current system. These improvements, 
based on exjoerience w ith the existing l.andsats, w-ill provide new capabilities in the space- 
craft, the sensor, the gramd system, and the overall system design. These system 
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capabilities - which emphasize improved vegetation analysis, prompt availability of data, 
freejuent coverage, and precise data registration and overlay for better elwnge detection 
will permit the Landsat D to capture already proven economic benefits in such diverse 
api)lications as: 

• Monitoring world-wide food prcxluctivity 

• Mapping agricultural land use 

• Monitoring rangelands 

• Surveying forest resources 

• Managing critical watersheds 

• IX-'tecting land use changes 

• t^l /mineral exjjloration 

An artist's concept of the Landsat I) system is shown in Figure 1-1- The spacecraft will be 
based on NASA's new Multi-mission Moduhir Spacecraft (MMS) and will operate two remote 
sensing instruments: a Thematic Mapper (TM), with 30 meter ground resolution, and a 
Multispectral Scanner (MSS), with 80 meter resolution. The system provides two data 
communication paths to the Earth; one is a direct readout link for ground stations Oxjth 
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domestic and foreijin) within range of the spacecraft, and the other is a relay link via the 
Tracking and Data lelay Satellite System (TDRSS) for nearly full glolxil coverage. The 
spacecraft will be in a sun-synchronous orbit with a descending node time of 9:;i0 AM 
(similar to current l.andsats). The orbital altitude and inclination will provide near global 
coverage of the land and near coastal regions with a repeat cycle every 16 to 18 days. 

The use of the new' MMS spacecraft r.s the basic bus will provide both improved sensor jx>int- 
ing accuracy (iO.Ol degree) and stability (10 degrees /second). These improvements will 
manifest themselves in more accurate and more straightfonvard geometric corrections of 
the image data; both relative (image to image) and absolute (with respect to the Earth's surface). 
The MMS incorjxjrates modular subsystems in the key areas of jx)wer, attitude control, and 
command and data handling. This modularity together with the comjxitibility for both conven- 
tional and Space Shuttle launches will enable in-orbit rejjair and refurbishment of the spacecraft. 

The Thematic Map|)cr, TM, is an evolutionary' improvement of the MSS and provides several 
significant capabilities. The spatial resolution on the ground has lx;en reduced to 30 meters 
(compared to 80 for the MSS) which will allow radiances to be measured for areas (pixels) 
less than one sixth the size as for the MSS. The TM will incorjwrate si.x spectral bands 
(and liave the capability for a seventh) which have Ix'en located primarily on the basis of 
their ability to discriminate vegetation (a fundamental application of remote sensing). In 
addition the radiometric sensitivity of the TM has been improved by reducing the sigTial-to- 
noise characteristic and increasing the levels of digital quantization. These sensor changes 
combine to cause the TM to have a data rate of 120 Mbps, (an order of magnitude increase 
over the 15 Mbps of the MSS). 

For remotely sensed multispectral data to-be truly practical tor many potential operational 
users (agricultural analysts, hydrologists, etc.) it must be received by them in usable 
form within 48 to 96 hours after imaging. Promptness in receiving data products is one of 
the most critical aspects of the landsat System. 

The l andsat D System will be thoroughly integrated with the needs of operational users. It 
will include improved preprocessing of all data, central data processing, archiving and retrieval. 


3 



low-cost receiving and data centers for large volume users (such as the U. S. Department 
of Agriculture) and provide maximum efficiency and economy in utilization by state, regional, 
and foreign users. Featuring the rapid electronic transmission of all data, the Landsat 
U system will reduce the time between satellite imaging and user reception of tlata 
to the required 48 to ‘.)G hours. 

As illustrated in the artist's concept the system provides tw'o data links to the ground. The 
first link, for both MSS and Thematic Mapper data, is directly from the satellite to domestic 
and foreign ground stations as the satellite jxisses through their reception areas. The 
second link is via the Tracking and Data Relay Satellite System (TDRSS). As shown, the 
data is transmitted to a TDRSS satellite, in stationary orbit, and relayed to the TDRSS 
receiving station. The TDRSS receiving station transmits the data via a domestic communi- 
cations satellite to a central data processing facility that, in turn, relays the data to any Icx-al 
data distribution center equipped to receive it. This link, via TDRSS and the communica- 
tions satellite, will thus have global acquisition and relay capabilities, providing rapid access 

to Thematic Mapper data for users throughout tlie world. Both data links have a planned 
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maximum data capability of 135 Mb/second it a 10 bit error rate. 


The Landsat D system described is currently in the planning stages by NASA. As part of 
the planning for this future system, NASA has undertaken a series of studies, with General 
Electric and others, to investigate various system options. This particular study is one 
of seven conducted by General Electric to exjilore different aspects of the total ground system 
that will Ije required by Landsat D in order to meet the overall mission objectives. 

'J he seven ground system studies are: 

1. Local User Terminal Study - an investigation into the requirements and 
options available for direct readout (primarily foreign ground stations) 
of Landsat D data. 

2. User Data Processing Study - an effort to estimate the sco)xj, size, and 

cost of the major user data processing system requirements. 
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3. Data Process in g; racility Study - a requirement and sizing study to provide 
preliminary estimiites of the scoix 2 and cost of NASA's central landsat D 
data processing center. 

4. GSIX Research & Development Study - a survey and analysis of the- functions 

and facility required of NASA to continue the basic research on spaceborne 
remote sensing and its applications. 

5. Operation Control Study - an analysis of the modifications necessary to 14 ) grade 

or modify the NASA Operations Control Center (OCC) for landsat D, 

6 . Data Transmission and Dissemination Study - an investigation into the options 
and limitations of various data communication alternatives including centraliza- 
tion versus decentralization. 

7. Position Determination and Correction Study - an analysis of the imixict and 
alternatives afforded by the MMS spacecraft of l.andsat D on image geometric 
correction. 


1.2 THE LANDSAT D GROUND SYSTEM 

A top-level functional diagram of the landsat D ground system is presented in Figure 1-2. 
The five major subsystems included are the Data Injxit Subsystem (DIS), the Central Data 
Processing Facility (CDPF), the Product Generation and Dissemination Facility (PGDF), 
the Data Management Subsystem (DMS), and the Agriculture Utilization Subsystem (AUS). 
Each of these subsystems is briefly described below. 
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Figure 1-2. Ground System Concept 
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The tkita Input Subsystem (DIS) receives 120-135 Mbps clata from the TDRSS via dedicated 
cable interconnection. The prime functions of the DIS are to record the raw infxit data, to 
perform cloud cover detection and scene editing, and to compute geometric correction 
matrices on a |X?r swath basis. 

The Central Data Processing Facility (CDPF) receives edited data from the DIS and performs 
standard operations to all data. These operations include radiometric correction and data 
reformatting to a band-interleaved-by-line (BIL) format. 

The Product Generation and Dissemination Facility (PGDF) is the maii; Interface between the 
I.andsat D ground system and general users. This facility provides l.andsat D data, in 
either digital tape or film format, to users on demand. The data, which may lx? geometrically 
corrected to various map projection systems or enhanced as requested by the user, is 
available in a variety of sizes, formats, and media. The PGDF also houses and manages 
the system archive. 

The Data Management tiibsystem (DMS) provides the central point of control and data base 
management for the I.andsat D ground system. Its prime functions include management of 
user demand, the system archive, system communications, and system redundancy. The 
DMS als(i maintains system status, production statistics, operations logs and administrative 
services. 

The Agriculture Utilization Subsystem (AUS) receives data directly from the CDPF and per- 
forms those operations necessary to produce world crop production forecasts on a periodic 
Ixisis. The operations to be ix?rformed include geometric correction, sample segment extrac- 
tion, multispectral analysis, and areal and statistical analyses. It is included here as part 
of the ground system because it represents the first major user of Landsat D dnta. 

Several other major subsystems included as part of the I.andsat D ground system were con- 
sidered. These include the Operations Control Center (OCC), the GSFC Research and 
Development Facility, and the Hydrologic Land-Use Utilization Subsystem. The OCC per- 
forms the functions required to plan, schedule, operate and evaluate spacecraft and payload 
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operations. The R&D Facility enables NASA to perform research related to the Uindsat 
program and its applications. The llydrologi c/Land Use Utilization Subsystem is similar 
in concept to the AUS and will generate land use maps over watershed areas within tiie US. 

1.3 POSITION DETERMINATION AND CORRECTION SUMMARY 

The l4indsat D system ^^^ll require a significantly more accurate method of determining the 
geometric correction functions for image data than is currently available . in additi(*i, 
more timely transmission of ephemeris data to the geometric correction process will lx* 
dictated by the operational nature of many of the anticipated user programs. 

Use of the Multi-mission Modular Spacecraft (MMS^ as the basic space vehicle for the 
Landsat D program will result in a considerable improvement in attitude control capability 
over current (Landsat 1 and 2) vehicles. In spite of this improvement, the vehicle attitude 
is still the overriding error source in the ability to determine the precise ground location 
of the sensor's instantaneous field of view. This means that predicted ephemeris, available 
24 hours ahead cf data acquisition may be used for the position of the spacecraft in all geo- 
metric correction calculations. 

The geometric accuracy of the Landsat-D imagery after simple corrections which model 
known error sources will be approximately 275 meters (RMS). Although this accuracy is 
significantly better than the present Landsat capability of 1-2 kilometers, the use of control 
ixjints to coiTect the imagery, on a swath basis, will be requ red to improve the {X)Sition 
and registration accuracies to approximately one-half pixel. 

In addition to improved accuracy, swath correction (possible because of very low* MMS 
attitude rates) results in about a one-hundred to one reduction in the number cf control 
|X)ints recjuired to be surveyed, stored in the data base and correlated. 

The rejxirt provides an assessment of accuracy of the knowledge of Landsat D siiacecraft 
ephemeris data, an evaluation of the impact of expected attitude and alignment accuracies 
and analysis of the various options for the combining of precision ephemeris and attitude 
data with scene image data. Finally, this rejxrt will characterize the potential geometric 
correction system in order to determine overall system costs and impact on other system 
e lements. 
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SECTION 2 


ERROR ANALYSIS 


2. 1 BASIC ASSl’MPTIONS 

Position accunicy is defined as the ability to relate a point in the image to its true position 
on the earth's surfaee as represented by some map projection. Registration accuracy is 
defined as the ability to overlay two images of the same ground point taken at different points 
in time. Throughout the following analysis, accuracies will be defined in terms of one sigma 
(or RMS) which refers to a 68.3 percent probability of occurrence. Since the positional 
effects of some error sources are pro|X)rtional to radial distance from the image center, RMS 
errors for such cases are represented by using the radial distance which delimits 68.3 jiercent 
of the image area. Wliere an error source causes both an X and Y positional shift, the vector 
sum is computed. 

The Multi-mission Modular Spacecraft (MMS) will maintain its attitude in orbit to within 
i0.01° (each axis) with a rate stability of ±10"^ deg/sec. The payload will consist of two 
sensors; a five-band Multispectral Scanner (MSS) and a six-band Thematic Mapper (TM). 

Each sensor has a 185 km swath width. The instantaneous field of view on the ground for 
MSS is 79 meters and for TM is 30 meters. 

2.2 ERROR SOURCES 

The effect upon position accuracy of each major category of error source is shown in Table 
2-1. This represents the system |X)sition accuracy after correction for earth rotation, initial 
alignment biases and sensor scan mcxieling. Compared to a 1-2 km error reported for landsat 
1 and 2 data, the Landsat-D system will perform three to seven times better even beforc control 
point analysis. 

2.2.1 ATTITUDE CONTROL 

The MMS is a very accurate and stable platform (see Reference !)• rhe RMS pointing specification 
(both knowledge and control) is tiO.Ol degrees for each axis, which is fifty times better than 
Landsats 1 and 2. Tlie RMS specification is 10”^ degrees /second for each axis, which 
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Table 2-1. Error Source Summary 


Error S(xirees 
Attitude Control 
Ephemeris 
Alignment 
CXher MSS 
TM 


Effect Upon Position Accuracy 
(After X & Y Travslction & Rotation) 

175. 2 Meters 

39.2 

206.4 

33.4 

24.0 • 


RSS Totals MSS 275.6 (3.5 pixels) 

TM 274.6 (9.2 pixels) 


is 15, 000 times belter tlian landsats 1 and 2. This performance will be acn>.'iplished by 
star tracker updates every 5-10 minutes to the inertial reference unit. 


Figure 2-1 shows the effect upon the image of each Attitude Control Subsystem, ACS, error 
source. The RSS total of all these error sources is l'^5.2 meters of position inaccuracy. 



Figure 2-1. Attitude Control Induced Errors 
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2.2.2 EPHLMEHIS 


An analysis was performe 1 by the Ojierational Orbit Support Branch at GSFC to determine 
the typical growth in radial, along track, and cross track prediction errors for Landsats 
1 and 2 (reference 2). The results were obtained by predicting an orbit obtained from a 
two day definitive arc into the best four cxinsecutive tw’o day definitive arcs and extracting 
the maximum compare differences ever days 1 through 8. The growth in prcMJiction error 
is a ftmetion of the inability to model the effects of both atmospheric variations and the 
variable number of attitude gatings. 

The RMS errors in predicted ephemeris after one day are showm in Figure 2-2 along with 
their effect upon ^x)sition of the image. Also shown are the RMS errors in best fit ephemeris 
taken from an analysis done for Landsat C (reference 3). Since the total system position 
accuracy will improve Iw’ only two meters if best fit ephemeris is used instead of predicted, 
the baseline Landsat D system will use predicted ephemeris. This approach permits real 
time correction of the data, while use of best fit ephemeris causes delays of at least two days. 
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Figure 2-2. Geometric Errors Induced by Ephemeris Uncertainties 
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2.2.3 ALIGNMENT 

The attitude control nuxiule v^ill be aliened to the MMS structure to within ±200 arc seconds 
(RMS) in each axis (see reference 1). It is reasonable to assume that the payload modulo 
will be similarly aligned to the MMS so that the initial alignment bias between ibe sensor optical 
axis and vehicle pointing vector will be ±200^2 st^. This bias will be measured before 
launch and can lx? removed from the image data. The alignments are expected to vary over 
the orbital {>eriod by as much as ±30\/2 s^ in each axis. This error source cannot be 
removed until control point analysis has mtxlelled its positional effect as a function of time. 
Figure 2-3 shows the alignment error sources by magnitude and their effect upon position 
accuracy. 
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Figure 2-3. Errors Intrrxluced by Alignment Uncertainties 
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2.2.4 OTHER ERROR S0URCP:S 

Figure 2-4 lists the remaining error sources which contribute small but non-negligiblc 
amounts of position accuracy to the imagery. Spacecraft velocity will Ix' known very 
accurately (reference 4), and its positional effect is not large. The knowledge of when the 
iniige was taken has been estimated very conservatively at ±1 msec. The vehicle altitude 
ccHild change by as much as 200 meters between ephemeris ujxlates causing a scale change 
in the cross-track direction as shown in the figure. The sensor scan is nonlinear and can 
be modelled so that residual errors are small. Internal source errors include such items 
as line start jitter, detector misalignments and scan-to-scan repeatability (see references 
5 and 6). 
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Figure 2-4. Errors IntrcKluced by Other System Uncertainties 
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SECTION 3 


CORRECTION PROCESSING 
3.1 IMPLEMENTATION SCHEME 

The errors present in the l^andsat D attitude, ephemcris and alignment data make it necessary 
to establish ground truth before more precise correction is possible. By determining the 
image coordinates of well defined geographic features to be used as Ground Control Points 
(GCP's) in the scene, the correspondence between image coordinates and the actual control 
point location on the earth's surface can be used to compute correction functions. 

Search areas will bo established throughout a swath where it is most likely that high contrast, 
cloud free, control |x)ints will be present in one Thematic Mapper band. The map coordinates 
of these search areas are adjusted for earth rotation, earth curvature, initial alignment biases, 
and sensor scan mtxleling. This is essentially an X and Y translation plus a rotation correction 
which reduces the location error to 275 meters. If a search area of 10 by 10 pixels is used, then 
the probability that the GCP is present within the search area will be 68.3%. 

The image coordinates of a ground control point are the sample and line numbex's in the digital 
data. There are many techniques which could be used to determine the image coordinates. 

Most automated techniques reejuire that subimages of the selected geographic features (window 
areas) have been identified in other Landsat scenes and entered into a ground control ixxint 
library’. One of the candidate techniques (e.g. Sequential Similarity Detection Algorithm 
(SSDA), Edge Detection or Fast Fourier Transform) will be applied to the image data in order 
to locate GCP's. SSDA employs a I'andom sampling sequence for comparing w indow area pixels 
with search area pixels; test results using this technique on Landsat data ;u:e available in 
reference 7. Fast Fourier Transform techniques compute the mathematical cross-correlation 
of the two images and register them where the correlation surface peaks (see reference 8). 

Edge detection and correlation are currently used to correct I.andsat images for the 1J\CIE 
project. 
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C)n the basis of the ixjsitions of GCP’s in the image, and their known geographic location, a 
least squares fit (or affine transformation) is comjxited, which results in a mapping function 
(global mapping ix)l\T»omial) to lx; used to com|xite the image coordinates for every pixel 
position in the corrected image. 

Since computation becomes excessive if the mapping |X)lynomial is computed for the more 
than 37 million pixels in one Thematic Mapper image, relatively few intenoolation grid points 
(anchor points) will be mapped with the full mapping function. The other out{xit image jx)ints 
will be located by bilinear interjxjlation from these anchor |X)ints. The interjx)lation grid is 
selected so that the errors in inter{xjlation are no more than 0.25 pixel. 

After a jxjint in the output plane has been mapped into the input image plane, a radiometric 
intensity must be assigned to it. In general, the distorted image plane will not have a radio- 
metric intensity corresponding to the jxjint mapped from the output plane. This requires that 
the input plane be resampled at the output point locations to obtain the correct intensities. 
Studies liave shown that cubic convolution resampling produces less error than either nearest 
neighbor or bilinear interpolation (see reference 9). 

Corrections may be applied to the image data as requested by the users, or data tapes may 
Ix' supplied with the correction coefficients as a header. The expected accuracies in the data 
after control point piaxessing are shown in Table 3-1, as a function of the number of control 
points correlated per swath. All errors are one-sigma (RMt^ values in terms of meters tn 

the ground. The GCP correlation error can vary from 0.3-0. 8 pixels depending u{XJn how- 
much the image is obscured clouds, haze or noise (see reference 8). At a USGS map 
scale of 1:24,000, National Map Accuracy Standards require that landmarks be displayed 
within 0.5 mm of their true position. The error in reading the true GCP coordinates fx’om 
a map is therefore 12 meters. Both of these errors will reduce by a factor related to the 
square root of the number of control points. 
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Table 3-1. Residual Geometric Errors After GCP Correlation 


Error Source 

Residual Errors (Meters) 


1 GCP 

2 GCP's 

4 GCP's 

7 GCP's 

Position Accuracy 





CtCP Correlation Error 

12.73-33.94 

9.0 - 24.0 

G.37 - 16. 97 

4.81-12.83 

Map Location of Control Point 
0. 5 mm at 1 :24000 

12 

8.49 

6 

4.54 

Residual Yaw Error 

13.4 

13.4 

12.19 

9.92 

Non-Linear Scan 

5 

5 

5 

5 

Resampling Error 
(cubic convolution) 

1 

1 

1 

1 

Grid lnter|X)lation 

4.3 

4.3 

4.3 

4.3 

Internal Sensor 

4.G 

4.G 

4.6 

4.6 

1 RSS 

23-39 

20-30 

17-23 

14-19 

1 Registration Accuracy 





i Remove GCP Map Location 

i 

12 

8.49 

G 

4.54 

1 RSS 

1 

1 

20-37 

18-29 

16-22 

13-18 


Simulations liave Ixsen performed to determine yaw by using GCP's (see Reference 10), and 
the results are reflected in Table 3-1 above. In that same reference the Thematic Maiper 
scan profile was assumed to be modelled to an accuracy of 0. 17 pixel (i.e. residual scan non- 
linearity' of 5 metei’s.) GE correction systems use cubic convolution (sinX)/X resampling 
with a resolution of typically 1/lG iiixel. The one-sigma resampling error is therefore 
1/lG \f3 pixels or one meter. The correction function for a scene is not computed for every 
pLxel but for an array of points uniformly spaced over the scene. The correction value for 
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any given pixel is obtained by a bilinear interpolation of the surrounding array points. 
ExiJcrience at GE indicates that iO.25 pixel is a realistic error having a uniform distribution 
(7.5/ \3 meters). In reference 6 the internal TM sensor errors such as start of scan stability, 
across scan non-linearity and detector misalignments combine to 4.6 meters. Wlien 
registering two images, the map coordinates of the control points are not relevant; thus 
the error due to map location is removed from the error budget shov n in Table 3-1 above. 


3.2 CORRECTION FUNCTIONS 

Certain error sources can be removed from the data since they are known or predictable 
without recourse to control point analysis. Removal of these errors results in the 275 
meter position accuracy described previously. The principal corrections are detailed in 
the following paragraphs. 

Images become skewed as a result of the rotation of the earth during the finite frame time. 

The linear velocity of the earth at format center is 

^ _ 2 radians/day x R x cos L 

86400 seconds/day 

where R = earth radius 

L = latitude of format center 

The inclination of the orbit to the local meridian at format center is 

I = sin‘^ sinE 
tosL^ 

where E = oi'bit inclination at eciuator 

The distortion components are therefore 

AX = V (cos I) At 
AY = V (sin I) At 

where At = time from format center 

Thus, earth rotation correction is predictable and involves simple X and Y translations of 
the data. 


16 


The sixicecraft axes are defined as follows; roll is coincident with the velocity vector, 
pitch is normal to the velocity vector parallel to the ground and yaw is the vertical axis of 
the spacecraft normal to the ground. Alignment biases in the roll axes result in the following 
image location errors: 

AX^ = H (tan - tan B) 

AY = 0 
r 

w'here 11 = S/C altitude 

B = angle between S/C nadir and instantaneous field of view 
0 = roll alignment bias 

The errors resulting from pitch alignment biases are as follow's: 

= 0 
P 

AY = H tan d 
P 

where d = pitch alignment bias 

Position errors due to yaw' alignment biases are; 

AXy = X(l-cosy) 

AY = X sin Y 

y 

w'here X = distance from nadir to the ground point in the image 
y = yaw alignment bias 

The sensor scan mechanism will deviate from the ideal linear case. Deviations from a 
lineiir scan profile can be determined empirically and incorporated into a scan model. 

The scan model can be applied to image data in order to reduce residual errors to a sub-pixel 
level. 

3..‘1 CONTROL POINT TRANSFORMATIONS 

As mentioned in [xiragraph 3.1a mapping function is derived fromthe positions of control 
|x)ints in the output and image plane. Two likely candidates for use as mapping functions 
are the affine transformation and the least squares transformation. 
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The affine transformation is a mapping which accounts for distortions due to translation, 
scale change, rotatitm, as|x>ct ratio and skew. To map a point from outjxit plane to input 
image plane requires three GCP's which form a triangle containing that point. The 
detailed algorithm is shown in Figure 3-1 (from Reference 11.) 
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Figure 3-1. Affine Transformation 
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The other method of mapping involves the use of least squares polynomials. These poly- 
nomials have the form: 

+C^X^^C 
6 

The coefficients for these coefficients are computed from the GCP data with the minimum 
number of GCP's increasing with the degree of the polynomials. Also the placement of 
GCP’s should be such that the comers and edges of the image area are covered to prevent 
large edge errors in the polynomial functions. 

3.4 MAP PROJECTIONS 

The user may be offered his choice of map projection, i.e. Universal Transverse Mercator 
(UTM), Space Oblique Mercator (SOM), and Lambert Conformal Conic (LCC). 

In a transverse Mercator projection the Earth is projected onto a cylinder whose axis falls 
within the equatorial plane. The line of tangency between the Earth and the cylinder is along 
a longitude meridian. The Universal Transverse Mercator takes the cylinder's axis and 
rotates it in six degree steps as shown in Figure 3-2 (Reference 11). Thus the cylinder fits 




Z = C ^C\X+C„Y+C„ 

o 1 2 3 4 5 
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Figure 3-2. Universal Transverse Mercator I^’ojcction 
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the Karth's surface closely over ^3 degrees from the longitude of tangency for each UTM zone. 
Surface coordinates are measured in meters from the central target meridian and meters from 
the equator. The Northing and Easting coordinates together \\lth the zone number define a 
point on the Earth's surface in UTM. Although it is not defined above 80° latitude, the UTM 
projection has the advtmtage of being a Cartesian coordinate system with identical units of 
measure along each axis. There is a disadvantage in terms of distortions due to projecting 
areas from one zone into the adjacent zone. This occurs at higher latitudes (above 60°) where 
the spacecraft heading angle increases at a rapid rate and the zone widths decrease causing a 
rapid traversing of zones. (For more detail see Reference 13). 

The Space Oblique Mercator (SOM) projection involves a cylinder whose circumference is 
tangent to the satellite ground track as shown in Figure 3-3. It is defined in Reference 12. 




Figui’c 3-3. Space Oblique Mercator Projection 
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The l^imbert Conformal Conic (LCC) is the projection of the Earth spheroid onto a cone whose 
axis coincides with the polar axis of the spheroid as shown in Figure 3-4. (from Reference 13). 
The cone is secant to the Earth, intersecting along two jiarallels of latitude which are on the 

same side of the equator. Meridians ap|x;ar as straight lines radiating from a |X)int Ijeyond 
the mapped area. l.atitude parallels appear as arcs of concentric circles which are centered 
at the point from which the meridians radiate. 



con««^<ir> pn*{«-rt«on 


Figure 3-4. Lambert Conformal Conic Projection 
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One slight problem with UTM and LCC projections is that the scan lines of image data will 
net be normal to the north vector. This is illustrated in Figure 3-5 for the case at the 
equator. In order to project the scan lines normal to a north vector it is necessary to rotate 
them through ever increasing angles as latitude increases. The cost of a rotational buffer 
increases rapidly with latitude also, as shown in Figure 3-6. The sa\1ngs in hardware cost 
outweigh the minor disadvantages of not rotating scan lines. Film products will simply be 
manually rotated to overlay on maps. Slight discontinuities in producing mosaics from high- 
latitude, non-parallel ground tracks will result. Digital map registration will require user 
applied rotations, but these should involve significantly less data and smaller buffer costs. 

Thus, the I.andsat D system will rotate scan lines for yaw and alignment biases only, which have a 
much smaller cost impact since the angle is approximately *0.10 degrees. 


SPACECRAFT 
GROUND TRACK 


EQUATOR 


NORTH 
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Figure 3-6. Rotational Buffer Cost 
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SECTION 4 


CONCLUSIONS 


Since the l^ndsat D Sj'stem will utilize the Multi-mission Modular Spacecraft, ACS pointing 
and rate errors are significantly less than those prcxluced by former Landsat spacecraft. 
The MMS stability permits control point processing and correction of the imagery on a per 
swath basis. Per swath correction minimizes processing complexity and the amount of 
hardware required. 

Since ephemeris errors are small compared to ACS and alignment inaccuracies and best 
fit ephemeris is not significantly better than predicted, the predicted ephemeris will be 
used. This offers advantages because correction processing is not delayed by' the best fit 
computations and can be available in near real time. 

A choice of map projections may be offered to the users since no one projection will satisfy 
all users. This approach is the most flexible way to insure user satisfaction. UTM and 
LCC map projections could require scan line rotations in order to orient them normal to a 
n<r th vector. However, the rapidly rising buffer costs incurred by large rotations is a 
strong deterrent. Sigiiificant cost savings coupled with a clear lack of universal user 
preference has driven the decision not to rotate scan lines for map projections. 
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